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Abstract. We propose a method to analyze the dynamics of systems exhibiting 
slow relaxation which is based on mesoscopic non-equilibrium thermodynamics. The 
method allows us to obtain kinetic equations of the Fokker-Planck type for the 
probability functional and their corresponding Langevin equations. Our results are 
compared with the ones obtained by other authors. 



1. Introduction 

The problem of the glass transition has attracted much attention in the past years, and 
has been subject of intense experimental, computational and theoretical investigations. 
Mode coupling theories and replica techniques have been able to capture some trends 
of the phenomenology of these systems. Although interesting theoretical results have 
been obtained by means of those theories, the understanding of the problem is far from 
being complete [l], ||]. 

Supercooled liquids and colloids near a glass transition are characterized by a 
extremely slow relaxation of the density variable. As the glass transition is essentially 
a dynamic transition, it becomes important to provide a dynamical description of 
the evolution of the density variable. In addition, the presence of slow dynamics in 
such systems emphasizes the importance of the activated processes, since they are the 
main mechanism able to make the system evolve and escape from its trapping in the 
metastable states. A complete description of these systems should then take into account 
such processes. 

In this paper we present a simple mesoscopic formalism able to derive kinetic 
equations of the Fokker-Planck type for the probability functional of the density fields 
in supercooled liquids and dense colloidal systems. This approach is based on the 
hypothesis of local equilibrium in the phase space of the system, which allows us to 
describe the processes leading to variations in the state of these systems by means 
of non-equilibrium thermodynamics. The extension of that theory to the mesoscopic 
level of description has been referred to as mesoscopic non-equilibrium thermodynamics 
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The plan of the paper is as follows. In Section 2, we derive the kinetic equation 
for the probability distribution functional of the momentum density and mass density 
fields. By adiabatic elimination of the momentum density we obtain the kinetic 
equation describing the slow dynamics, giving the configurational changes in the system. 
In Section 3, we introduce fluctuations of the distribution functionals obtaining the 
corresponding Langevin equation. Finally, in Section 4, we conclude summarizing our 
main results. 



2. Kinetic equations 

We consider the system, either a pure liquid or a colloidal suspension, as a continuum 
divided in tiny cells of volume Vq. In each cell we define a mass density p(r) and a 
momentum density g(r) = pv, with r being the position vector of a cell and v the 
velocity field. For simplicity, we will introduce the compact notation {£} = {p, g}. 

Let us consider P({T} ,t) as the distribution functional in phase space, normalized 
according to the relation 

P({T},t)6F = N, (1) 

where N is the number of possible states in phase space. Related to this distribution 
functional we introduce the phase space entropy S(P({T} ,t)), which satisfies the Gibbs 



equation introduced by mesoscopic non-equilibrium thermodynamics fllO 



AS= -\ //*({T},t) A P({T},t)6T. (2) 

Here T is the temperature, and p represents a chemical potential defined in phase space 
satisfying 

SS 

fi=^. 3 
5P 

Its expression can be obtained with the help of the Gibbs entropy postulate [|I~0|1 

S = -k [ Pin -J^—5T + S Leq \ (4) 

J pl-eq. ~ V ' 

where k is the Boltzmann constant. In this expression 

pi rn- = exp { _p [/io + Hk {E } + Hu {p}]} (5) 

is the distribution functional at local equilibrium, with (3 = 1/kT, and p Q the chemical 
potential at local equilibrium. Moreover, we have defined the kinetic energy functional 

and the potential energy functional 

Hu {P}=jf {P} dr (7) 
where / {p} is the free energy density functional. 
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By taking variations in eq. ( |) we can obtain 

A S = -k [ APln -J^—5T + AS Leq -, (8) 

J pl.eq. — v ; 

where 

A S Leq - = A N = -~ J fx A P5T. (9) 

Thus, after comparison of eqs. (§) and (|j) we can infer the expression of the chemical 
potential 

P 

fi = fi + kT\nj^. (10) 

We now assume that the distribution functional evolves according to the continuity 
equation 

where the dot over the field variables indicates time derivative and J g is a diffusive 
current in momentum space. 

>From eqs. @ and (|Tl| ) we can obtain the rate of the entropy variation that is 
written as 

In this expression, we can identify the first term on the right hand side as the rate 
at which entropy is supplied to the system by its surroundings through the external 
constrains, given by 

deS 

dt T J J n \5p r ~ ' 5g 

The second contribution corresponds to the entropy produced inside the system due to 
the irreversible processes, whose value is 

E = -?// J .-5 <Wr " (14) 

where a partial integration has been performed. This entropy production has the usual 
form of flux-force pairs from which we can infer the phenomenological relation 

J 9 = -yUr,r^)'^, (15) 

where the Onsager coefficients L satisfy the Onsager relations L(r,r') = L{r' ,r)\ in 
which t stands for the Hermitian conjugated. By computing the functional derivative 
and assuming locality in the coordinates, i.e. L(r, r') = L(r)5(r — r'), one obtains 



- / U[-pP+ — .gP)drSr. (13) 
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where H_{r) = L{r)/T P can be interpreted as a mobility tensor. Its expression follows 
from the Navier-Stokes equation [[□]] 

H{r) =- v (^VV + IV 2 ) " f V V (17) 

with i] and £ being the shear and bulk viscosities, respectively, and 1 the unit tensor. 

By introducing the current we have obtained in eq. (|i~6|) into the continuity equation 
(fLl|) and using the expressions of p and g given by the reversible part of the Navier-Stokes 
equation |TT|, this yields 



dt J \5p^ 9 5g 



5Hu gg 

dp p 

'g-1 8 | 6H K 
$g Sg 



Pdr, 



(18) 



which constitutes the functional Fokker-Planck equation for the probability distribution 
functional P({F} ,t). 

In the diffusive regime, when equilibration in momentum has occurred and only 
configurational changes are possible, this equation simplifies considerably. To achieve 
the corresponding equation we define the reduced distribution functional 



and the current 



P5g 



gP5g. 



By applying the time derivative to eq. fll9|) and using eq. fll8|) we obtain 



OP 



dt 
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= / T V -Jpdr. 



S P 



(19) 



(20) 



(21) 



The probability current defined through eq. (|20|) evolves according to 



dJj, 

dt 



f dP x 



5P 



5p(r 



^g( r ) V -g(r')5gdr 



PV^Sg 

P 



pPV 



6 H } 



u 



5p 



-H_ ■ J p , 

P 



(22) 



where we have used eq. (|18|), and partial integration have been performed taking into 



account that Sg(r)/8g(r') = 5{r — r'). We now assume that -H_ ■ J ' p ~ r 1 J p , where 
t is the characteristic time scale of the inertial regime. In the diffusive regime when 

(23) 



t <C t, assuming equilibration in momentum we can write |12 

P = ${ r}P({ p },t), 

where 



$ = Z {pY l exp ^-/3/2 J ^~dr 



(24) 



Mesoscopic approach to the slow dynamics 

(9) P , 
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being (g) p = p / gP5g the conditional 



(— An 


~p{r) 


dr\ 


\2v J 


. Po . 





(25) 



with p being the uniform mass density of the system. Thus, by substituting eq. (|23|) in 
eq. (^), performing momentum integration and eliminating inertial terms we achieve 





with 



H {p} = Hu - kT In Z = E v - 



3kT 
2^7 



j In 


'p(r)" 







where we have used the result / g(r)g(r')Q5g = kTlp(r)5(r — r') 
substitution of the current J p , given in eq. ( |26|) into eq. ([H]) we obtain 

5 



(26) 

dr-kT In Z , (27) 
Finally, by 



<9P 



-fcT 



<5p(r-; 



[v-pWv] 
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<5p(r) 5p(r 



(28) 



which constitutes the functional Fokker-Planck equation in the diffusive regime, in 
which rkT plays the role of a diffusion coefficient. The stationary solution is the 
Boltzmann distribution functional P st ~ exp (— (3H) thus satisfying detailed 

balance. This equation has been previously obtained in [0 by means of projection 
operator techniques. It should be mentioned, as already been pointed out that the 
presence of a diffusive term in this equation incorporates the existence of hopping or 
activated processes. 



3. Fluctuating kinetic equations 

Our purpose in this section is to study the dynamics of the fluctuations in the 
distribution functional around a given reference distribution. To this purpose, we will 
consider that the distribution solution of eq. (|1|) corresponds to the average value over 
an initial distribution in phase space. Consequently, the actual value of the distribution 
functional will differ from the solution of eq. fll8|) in the presence of fluctuations. 

In our derivation of the Fokker-Planck equation, we will incorporate such 
fluctuations by applying fluctuating hydrodynamics in phase space |o|. Thus, we will 
split up the current J g into systematic and random contributions 

j, = jg + (29) 



The former is given precisely by the linear law ( pT5| ) whereas the latter defines a Gaussian 
white noise stochastic process of zero mean and fluctuation-dissipation theorem given 
by 

( J* ({£} , r, t)J*({T!} , r', f )>£ = 2kL5({L} - {r'})<5(r - r')5(t -t), (30) 



Mesoscopic approach to the slow dynamics 



6 



where L has been defined in eq. (pT5|), and P is an initial probability distribution. 

When employing the decomposition (|29|) in the continuity equation ([TT|), we obtain 

where P* is the fluctuating probability distribution. This equation can also be expressed 
as 



dP* 



5p 



V-9 



8 



Sg 



P V 



SHv 
bp 



V 



r 1 - + ^ 

8g 5g 



OR. 

P 



P*dr + I —■ J*dr, 



5_ 
Sg 



(32) 



which corresponds to the fluctuating functional Fokker-Planck equation. 

Using the same adiabatic elimination procedure performed in the previous section 
we can obtain the fluctuating functional Fokker-Planck equation in the diffusive regime, 
yielding 



dP* 



-rkT 



S 



8p{r) 



>From this equation it is possible to get the Langevin equation for the first moment 
of the distribution functional, defined through 

(34) 



P*Sp{r)+r,{r,t), (35) 



p(r,t)= J p(r)P*({p},t)6p(r). 

By multiplying eq. (J33|) by p{r) and integrating one obtains 
dp{r,t) 



ot 



rkT J [ V ■ p(r)y] 



5 



5H 



5p(r) Sp(r 



where r/(r, t) = — r y ■ J J^5p5g is the noise term, whose correlation, after eq. (|30| ) is 
used, is given by 

(rj(r, t)ij(r', t')) = 2kTr V r -p{r, t) V w S(r - r')5(t - t). (36) 



The first term inside the bracket in eq. 



is the noise-induced drift 14 , which arises 



from the nonlinear character of the dynamics which is due to the density dependence of 
the kinetic coefficient. 

Assuming P* = 5 ({p} — {p}), eq. (p5|) yields 



dp(r,t) , , 6H({p(r,t)\) . , 



(37) 



9t v rv , . v Sp(r,t) 

which coincides with the stochastic equation in terms of the number density field p(r,t) 
in Refs. [0, |15|]-[19|]. The average over the noise of eq. ([37]) yields the mean field 
equation proposed in Ref. |2( 
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4. Conclusions 



In this paper, we have proposed a method to obtain Fokker-Planck equations for the 
probability distribution functional of the density fields for a liquid in the framework 
of mesoscopic non-equilibrium thermodynamics fl7j. Keeping the essentials of non- 
equilibrium thermodynamics, and extending its range of validity to the mesoscopic 
domain we have proposed a Gibbs equation where the entropy depends on the probability 
distribution functional in phase space. The peculiar dynamics of glassy systems enables 
one to assume the existence of two dynamical regimes, one related to the fast processes 
or inertial regime, the other related to the slow relaxation processes or diffusive regime. 
After, adiabatic elimination of the fast degrees of freedom 0, we have obtained the 
Fokker-Planck equation for the probability distribution functional of the fluid mass 
density, which describes the slow dynamics of the system. The previous description is 
independent of the underlying microscopic model. The important requirement is the 
existence of well separated time and length scales in the system, which sustain the 
validity of the local equilibrium hypothesis in phase space. Supercooled liquids and 
dense colloidal suspensions are known to satisfy this requirement. 

At mesoscopic level we have considered the hydrodynamic fields of the liquid 
as stochastic variables whose deterministic dynamic is governed by the Navier-Stokes 
equations. From these equations we can identify the phenomenological coefficients in 
the kinetic equations. We have also analyzed fluctuations of the distribution functional 
by adding a stochastic contribution to the current in phase space which satisfies a 
fluctuation-dissipation theorem. From this procedure we obtain a Langevin equation 
for the mass density. 

One of the advantages of the formalism we have proposed in this paper is that the 
dynamic description in terms of a probability functional provides a complete information 
of the slow dynamics of the system. In addition, its validity is not restricted to the 
particular situations addressed in this paper. The same kind of description is perfectly 
applicable to a wide variety of different situations, among which we can mention the 
cases of granular materials |2lJ or nucleation and phase transitions in inhomogeneous 
media l|23|, 
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